Method of obtaining an anticipatory estimate of a cell&#39;s wireless coverage

ABSTRACT

The invention concerns a method anticipatorily estimating the radio coverage by a base station of a cell ( 1 ) of a cellular wireless telephone network, whereby  
     a numerical map ( 31, 32, 33 ) is devised in the form of a direct-access data matrix of local and independent specifications of the positions and of kinds of salients ( 41, 42 ) of a plurality of predetermined meshes ( 34 ) of a mesh topology corresponding to the map ( 31 ), and  
     whereby, in an operational stage,  
     a sampling beam rendering the radioelectric propagation conditions is fictitiously transmitted in the cell ( 1 ) along an initial path segment ( 51 ) from an initial position and along a determined direction and under determined conditions of propagation,  
     the position of the instantaneous site is compared with the mesh topology to identify an incidence mesh at a salient ( 41, 42 ), and  
     the post-incidence propagation conditions are computed on the basis of the data of local and independent data specifications relating to salients of the incidence mesh.

[0001] The present invention relates to a phase defining cellular, wireless telephone networks to be set up in a territory and in particular to a method of obtaining an anticipatory estimate of the territory's radio coverage for the purpose of determining operational parameters and the optimal positions of base stations or repeaters of the network, namely the boundaries of the corresponding cells.

[0002] As is known, a wireless telephone network includes a plurality of ground base stations which are interconnected by a wired telephone network and which can be accessed by mobile terminals when the latter are within the station's wireless cell.

[0003] Radio propagation within a cell must meet two essential requirements, namely transmission at less than excessive power by the base station and reception of sufficiently powerful signals at the terminals.

[0004] In the first place the range of each base station must be adequate to extend into an adjacent cell in order to preclude any danger associated with a loss of contact that occurs when a mobile terminal moves into another cell. Accordingly the power at the transmitter must exceed a rigorous minimum.

[0005] In the second place, because the radio links propagate substantially straight at ground level, the station's shadow zones caused by the local topography or by buildings must also be covered. A shadow zone is a zone wherein tile attenuation of radio propagation between a mobile terminal therein and a station falls below a specified sensitivity for radio circuits, whereby a received level is inadequate to properly detect transmitted bit packets representing voice or data to be exchanged. On the other hand, tile levels of transmitted power cannot be increased.

[0006] The reason is that, at the base station, any increase in power also would increase cell size, entailing undue interference between adjacent cells. Maximum power at the mobile terminals is limited on safety grounds and the operating time of charged batteries at the mobile terminals.

[0007] Furthermore tile number of base stations or repeaters must not be needlessly multiplied to take care of the microcells—the shaded areas—so costs and interferences can be managed.

[0008] In the prior art, the attenuations are ascertained at a plurality of sites in a cell by using a vector database from the National Geographic Institute (IGN) for instance, which represents a map of the geographic zone under consideration containing buildings and other above-ground structures. Different code words define the kind of above-ground structures, for instance woods, lodgements, water, which are specified in Lambert coordinates and in height above the local ground altitude relative to the sea.

[0009] The transmission of radio signals is simulated to provide an anticipatory estimate of the attenuation at any point within the cell. Propagation is modeled by computing a set of radio coverages in a computer as an electromagnetic beam centered on the base station, said propagation being in a given direction, and the cell propagation conditions are then calculated within the element of a solid angle subtended by the beam. Except for the propagation in free space when the station is in direct line of view with a simulated radio terminal—from which follows a known linear propagation attenuation—the beam path strikes obstacles which attenuate it or additionally deflect it, in particular within micro-cells where the stations are situated at an altitude lower than the buildings' roofs.

[0010] Illustratively, in a street, the beam can deviate by reflection or refraction. As a result tie aperture of its solid angle can increase.

[0011] Such calculations are repeated for a plurality of elementary solid angles distributed within a global solid angle of view of the entire cell, for instance a substantially horizontal annulus, in order to sample the various propagation conditions of the cell's space.

[0012] At each point of each beam's path, an operator consults tile vector database to determine whether an obstacle is present. The related vectorial calculations require considerable computing power and might entail a day's work, and in practice they must be initiated during the evening when conventional computers are used.

[0013] The objective of the present invention is to reduce the computing power required to calculate the radio coverage of such cells, whether the cells be large, i.e macro-cells, or microcells.

[0014] For that purpose the present invention relates to an anticipatory estimation method applied to tile radio coverage of a cell of a wireless telephone network, said estimation being carried out by a cell traffic-managing radio station using a topographical map as the database containing the positions and the nature of the salients of tile cell, and in this method a computing system, when in its operational mode:

[0015] simulates irradiation of the cell with a sampling beam representing tile radio propagation conditions along an initial path segment from an initial position and along a given direction and under specified propagation conditions,

[0016] by means of the computer reading the database, compares tile path segment to the data of tile topographical map in order to identify the position and the nature of any instantaneous incidence site on the path segment comprising a salient,

[0017] determines new conditions of propagation on a downstream path segment beyond said incidence site depending on the information about tile particular salient stored in the data base,

[0018] iterates the two previous steps a given number of times at other downstream path segments as called for,

[0019] determines the cumulative attenuation at any point selected along said point in light of all the conditions of propagation along tile full path, and

[0020] repeats all the above stages a plurality of times for a plurality of initial directions in order to sample the entire cell and in this manner determine an attenuation map of the selected sites,

[0021] where this method is characterized in that, in an initial phase,

[0022] The computing system devises tile database in tile form of a direct-access data matrix containing the local and independent specification of the positions and natures of the particular salients of a plurality of predetermined meshes of a mesh topology corresponding to the map, and

[0023] The system stores the geographic-orientation data of the salients in a matrix database (32, 33)

[0024] During an operational phase, the computing system:

[0025] Compares the position of the instantaneous site with the mesh to identify a mesh containing an incident beam,

[0026] Computes post-incidence propagation conditions from the specified local and independent data of the salient of the incident mesh, and

[0027] Computes a direction of a beam reflected from the downstream segment (52) based on orientation data.

[0028] The propagation data have been divided into a plurality of independent data, and hence have reduced size, they can be quickly accessed to be read in useful form.

[0029] Accordingly computation of the propagation conditions downstream of the impact site is based solely on local data and on direct access rather than sequential access as in the case of a vector database, said data defining, as regards the incidence site, the anisotropic radio propagations that determine, by attenuation, reflection or diffraction, possibly with scattering, and a possible angular deflection and new propagation attenuation.

[0030] In the absence of an obstacle including a salient, the conditions of propagation in free space are very well known. In particular in this case, tie data readout from a crossed mesh at once show there is no obstacle and therefore the next mesh is considered without necessity of the prior art's ponderous computations of reconstituting a local data vector based on a global database.

[0031] It must be borne in mind that herein the term “salient” in this document always means a propagation obstacle, including even a horizontal ground, that might at least partly absorb or reflect the radio beams.

[0032] Advantageously the computing system in its initial phase stores the geographic orientation data of the salient in a matrix database and, in operation, and, based on the orientation data, it computes the direction of a beam reflected from the downstream path segment.

[0033] Preferably for this case, storage is restricted to azimuthal data and, during operation, the direction of the beam reflected from the downstream path segment is computed by assuming that the reflecting salients are vertical.

[0034] As a result computations are limited.

[0035] In particular, data specifying the kind of above-ground structures are integrated into the matrix database, and, in operation, the computing system computes the propagation conditions on the path segment as a function of the type of above-ground structure.

[0036] In this manner, higher grade propagation computation is attained.

[0037] In order to take into account details of the salients, the data concerning the edges of the salients can be initially integrated into the matrix data base and, during operation, tie computing system computes a downstream-segment direction of the refracted beam from the edges-data.

[0038] Preferably again, attenuations from the salients are initially stored in the matrix database and, in operation, the computing system computes a propagation attenuation of the downstream path segment based on said attenuation data.

[0039] Again, in a similar case, the attenuation data can relate to reflection from the salients and be used to compute the attenuation of the downstream path segment and/or the attenuation data concerns propagation through the salients and are used to compute the attenuation of the beams passing through them.

[0040] Advantageously the propagation-attenuation data through the salients also include transition data between propagation media that specify salient-penetration attenuations in these salients, said latter attenuations are used to determine local penetration depth.

[0041] In order to attain higher quality modeling of the propagation conditions, an algorithm is provided to compute the beam's angular dispersion following incidence, tile computing system computes a plurality of directions of downstream path segments at specified attenuations that constitute a solid dispersion angle of the beam beyond incidence.

[0042] Also the consecutive incidences along the beam path may be counted and, as regards the second incidence, in order to compute tile farther out propagation conditions, the computing system assumes the beam is polarized at the first impact point.

[0043] In order to keep track of the magnitudes of the required computations, the computing system counts the successive incidences along the beam path and compares the total to a high threshold value in order to stop performing in the operational phase when the threshold has been reached. Alternatively, the computer system determines at each incidence point the cumulative attenuation and compares the latter to a maximum attenuation threshold value to cease carrying out tile procedural steps when the threshold is reached.

[0044] To optimally make use of the results attained, the attenuation map is stored in 3D.

[0045] In this manner the anticipated quality of the radio links according to the building floors can be estimated.

[0046] In particular, in order to develop the database, the computing system represents the map by a bundle of vertically extending pixel strings which it divides to forth stacked elementary volumes of meshes each including its particular data.

[0047] Preferably the computing system assumes that the initial position is at the station, though it also can assume the initial position is an arbitrary one within the cell, and the beam transmission direction is selected as a function positions and types of the nearby salients in order that the beam travels near the station.

[0048] After the propagation conditions within a microcell in contact with a cell—where tile latter is larger and also is termed “macrocell”—have been computed in the manner of the present invention, the computing system can compute the propagation conditions within the cell, and then smooth the two computations relating to a boundary zone between the cell and the microcell.

[0049] In this manner the method of the invention allows improved functional integration of the two kinds of cells.

[0050] The invention is described in the description below of a preferred embodiment mode and in relation to the attached drawings:

[0051]FIG. 1 is a geographic relief map constituting a database onto which is entered tile position of a cell of a cellular wireless telephone network being set up,

[0052]FIG. 2 is a radio attenuation function within the cell along a beam path including obstacles, the base station forming one end of the said path,

[0053]FIG. 3 is a top view of a building contour of part of the map,

[0054]FIG. 4 is similar to FIG. 3, including reflection and refraction at a building, and

[0055]FIG. 5 is a vertical section including propagation between buildings with in one macrocell.

[0056]FIG. 1 shows the position of a cell 1 of a cellular wireless telephone network designed for terminals such as denoted by 21, said cell position being indicated in a portion of the geographic map 31 also showing the intended position of a base radio station 11 managing the cell 1. The map 31 enables a determination to be made of the expected contours of a plurality of cells constituting a radio network to be set up by adjusting the number of cells, their sizes and position in order to optimize the equipment bulk while assuring the desired radio coverage at a specified service quality.

[0057] Actually the reference 31 in this instance denotes elements of the terrain, on and above ground, in the zone of interest. The corresponding geographic data enabling this display to be generated are stored in a salients database 32, 33 of a computer 30. This database comprises a memory module 32 specifying the shape and topography of the ground and above-ground terrain and is associated with a memory module 33 storing the ground and above-ground morphology and specifying the nature or features of radio propagation of the various localized ground functions within a given frequency range corresponding to the frequencies used by the base stations, said morphology data of the memory module 33 representing the various sites of the terrain under consideration together with the salients or shapes of the geographic memory module 32. Be it borne in mind that the term “salient” is construed broadly to denote any obstacle on which tile beam from the station 11 is incident directly or after having been deflected. Accordingly the term globally deals with the ground and above ground features. Besides the buildings, hills and the like, it also can denote the slopes of valleys or plains or expanses of water.

[0058] As is discussed further below, the data from the memory module 33 enable determining—directly or using a correspondence table involving the nature and data of radio propagation features—the perturbation imparted to an incident radio beam in order to ascertain the direction and the attenuation of tile amplitude of a corresponding downstream beam. In this instance the map 31 per se only serves a didactic purpose because the data defining it are contained in the computer-processed memory modules 32, 33.

[0059] In an alternate embodiment, the base station 11 is replaced by a station operating in the same manner but having a shorter range for the purpose of defining a microcell. As initially indicated above, the microcells are in zones of strong salients or in urban areas to cover the radio shadow zones of conventional cells.

[0060] It is assumed at this point that a narrow radio beam linking the station 11 to the mobile terminal 21 within the cell 1 encounters obstacles 41 and 42, respectively a building and forest trees. For simplicity of exposition, it is also assumed that the obstacles 41,42 of FIG. I do not deflect the beam path which therefore remains straight and undergoes no reflection or refraction by the obstacles.

[0061]FIG. 2 is a plot, with the ordinate in dB of the radio-signal level S as a function of the distance X covered along path of the abscissa. Starting at the transmission site, the signal level decreases with distance of propagation and therefore represents attenuation. This attenuation is the sum of the attenuations of various path segments, each of which corresponds to a specific propagation medium.

[0062] In this instance there are five segments consecutively referenced from 51 through 55 and respectively corresponding, as regards the first segment 51, to the attenuation along the air path from the base station 11 to a building 41; as regards the second segment 52, to the attenuation due to passing through tile building 41; as regards the third segment 53, to the attenuation due to the air path to the edge of a wooded area 42; as regards the fourth segment 54 to the attenuation due to crossing the wooded area 42 and as regards the fifth segment to the attenuation due to the air path to tile mobile terminal 21.

[0063] The attenuation per unit length of the radio propagation in the free space of air is shown by tile slope of the signal level on FIG. 2 and it is a well known physical constant at a given carrier frequency, as a result of which—and based on tile propagation distances contained in the database representing tile map 31—the computer 30 therefore is able to compute the three respective attenuations. On the other hand, the obstacles 41, 42 provide increased attenuations. Moreover, and as already mentioned, some types of obstacles also can deflect the path as shown in FIG. 4. To verify the anticipated radio coverage, tile field as it spreads from diverse sites through the entire cell 1 must be computationally estimated for the links between the base station 1 and the movable terminal 21.

[0064] The anticipated link budget so computed—i.e. the sum of the attenuations relating to the path segments 51 through 55—cannot exceed the difference between a maximum transmission level Nm of the base station 11 and a predetermined sensitivity threshold Sm of the terminal 21. This constraint also applies to tile up-direction of communication from the mobile terminal 21 toward the base station 11, in this instance 2 watts, and a sensitivity level of the station 11. These sensitivity levels take into account a detection code for propagation errors and for a bounded number of error bits in the bit packets exchanged through time-division channels of a radio frame.

[0065] For efficiency, the computer attenuation computations call only take a limited time, much less than one day's work.

[0066] For that purpose, tile map 31 assumes the form of a geographic relief map representing the cell 1 and the database 32, 33, and it further specifies tile positions of salients, such as 41 and 42 and others which intrinsically affect radio propagation conditions, such salients for instance being high buildings, woods, lodgements, lakes and others. Predetermined meshes are assigned to the map 31 in an initial stage, and the numeric data of the salients which the map comprises are assigned to each mesh 34 of tile geographic matrix thusly set up in order to make available a matrix database, i.e. a mosaic, of salient specifications, which can be used in tile operational stage.

[0067] In practice, a human operator or the computer 30 defines the mesh topology and loads into the memory modules 32 and 33 respectively the shapes of the salients of each mesh 34 and tile data describing the nature of such salients. As mentioned above, these data identify the salients such as woods, buildings, and a general correspondence table allows determining therefrom tile propagation conditions of the beams incident on the particular salient. As a variation, the propagation-parameter values are loaded directly into the memory module 33, without need to store the kinds of salients.

[0068] To anticipate an estimate of the radio coverage of cell 1 by the radio base station 11 managing the traffic of this cell 1, computer 30, operates the database 32, 33 to form tile relief map that specifies the positions and kinds of salients, such as 41 and 42, of cell 1.

[0069] The computer simulates transmission in the cell of a sampling beam representing the conditions of radio-electric propagation along an initial beam path, for instance 11, 41, from an initial to position and in a direction and propagation conditions that are fixed, that is in practice through air.

[0070] By reading the database, the computer compares the path segment 11, 41 with the relief map data, in this instance specifically those of the memory module 32, in order to identify the position and nature of any instantaneous incidence site in the path segment comprising for instance a salient 41.

[0071] Depending on the specific nature of the salient 41 under consideration, the computer ascertains new propagation conditions along a downstream portion of the path 41, 42 (segment 52) beyond tile incidence site at 41.

[0072] Where called for, the computer iterates tile two preceding steps a given number of times for other incidences at the downstream path portions 53, 54.

[0073] Depending on the propagation conditions over the full path from 11 to 21, the computer determines a cumulative attenuation at any point along tile path and it repeats all tile above steps a plurality of times for a plurality of initial directions in order to sample tie full cell and in this manner to establish an attenuation map of the selected sites.

[0074] Moreover—and having previously established during the initial steps the database 32, 33 in the form of a matrix or mosaic of data relating to tile local and independent specification of the plurality of predetermined meshes 34 of tile mesh topology corresponding to the map 31—it is possible therefore, as regards the operational stage:

[0075] to compare the position of the instantaneous site to the mesh topology in order to identify-the incidence mesh 34, and

[0076] to compute the propagation conditions after incidence (for instance the path segment 41-42) from the data of local and independent specification of the salient of the incidence mesh 34.

[0077] In general therefore the database 32, 33 specifies within each mesh, in particular in the memory 33, the local radio anisotropy, that is the distortions imparted to the radio beam such as deflection, attenuation, diffraction, polarization and others. In a certain way, for each mesh 34, an ellipsoid of anisotropy that determines the propagation conditions in space, is thereby defined, this ellipsoid relating to three inherent directions, for instance three orthogonal unit vectors having an abscissa and ordinate, such as local parallel, meridian and vertical.

[0078] Actually there are multiple ellipsoids as discussed because specifying the values of several propagation variables, makes it possible, for instance, to compute the exit direction of an incident radio beam as functional of its angle of incidence on the map 31, or else the attenuation corresponding to crossing the mesh 34 is a function of two directions, namely angles of incidence and exit. Accordingly this is a matrix that transforms the propagation conditions of the diverse meshes where the data are stored in independent zones of the memory module 33.

[0079] When the appropriate zone for the mesh 34 is read from the memory module 33, it is possible to rapidly ascertain the propagation conditions at the mesh exit point from the path segment at the entry point, for instance path segment 11, 41. In particular, one bit per mesh 34 can specify whether the mesh 34 under consideration does or does not contain an obstacle. By direct readout of the obstacle, which where called for is arrayed with its homologs in a rapid-access and compact register, absence of obstacle is detected immediately and the computer 30 at once moves on to examine the following mesh without computing new propagation conditions. Accordingly except for reading the bit indicating that there is an obstacle, the computer 30 then does not consult the data module 33 that specifies the nature of the radio obstacles. Calculation of the attenuation at the instantaneous site from mesh to mesh along the beam path is not required: The cumulative attenuation must be computed only when there is an obstacle by computing tie distance between tie two end meshes 34 of the beam path under consideration. Experience has shown that the computations applied to a 500 m mean-radius cell require a time of about 1 minute for a computer having typical computing power.

[0080] Having fixed or computed the angle of incidence of the beam in the space of the map 31, the computer 30 thereby can directly read the memory module 33 for all corresponding values of the propagation parameters of the beam exiting the mesh 34—provided the beam in fact is able to exit the mesh.

[0081] Considering that the map 31 is a relief map, the mesh topology of the data preferably and, as in this example, is implemented in three independent dimensions like those cited above. In other words, a bidirectional topology of so-called horizontal meshes can be defined, each horizontal mesh 34 being associated with a vertical extension volume divided at various altitudes that can be specific to each horizontal mesh 34 by means of planes or other surfaces in order to define volume elements, each storing particular propagation data in a zone of the memory module 33. In most cases, two volume elements suffice for one horizontal mesh, tile one below illustratively containing all of a building and the one above corresponding to free space. On the other hand as regards overhanging salients, such as arched buildings or bridges, a freely propagating third volume element must be provided underneath a volume element containing the obstacle under-consideration.

[0082] In other words, the map 31 is shown as a bundle of pixel strings for each mesh 24 and the vertically extending pixel strings are divided to constitute stacked volume elements of 3D meshes, each having particular data.

[0083] The specification data of the features or nature of the salients of each mesh 34, which were integrated during tile first stage, can correspond to one or more of the following data.

[0084] The specification data of tile salients in the memory 32 can include tile geographic orientation data of the salients, illustratively indicating a radio reflecting plane. Knowing the incident path portion or segment, the direction of the downstream reflected path segment is then computationally inferred said reflected segment being symmetrical relative to the normal at the reflecting plane and the site of incidence. In air, in the absence of an obstacle in the first Fresnel ellipsoid (direct propagation), the attenuation in the near field is about 20 dB/km over the first 500 m of beam path; beyond this attenuation rises to 30 dB/km.

[0085] Accordingly FIG. 3 is an illustration of the contents of a table of matrix data of tile memory module 33, the contents being schematically shown in graphic form in a part of map 31 for better clarity of exposition.

[0086]FIG. 3 includes several horizontal meshes with plotted building surfaces. The particular building under consideration has four straight walls 61 through 64 and a trapezoidal shape with four wall corners or edges 65 through 68.

[0087] In the memory module 33 the four meshes including one of the wall edges comprise data specifying that characteristic. Moreover the value of tie angle of the edge and even the orientation of its sides can be specified. Tile other meshes 34 crossed by one of the walls 61 through 64 include data specifying this feature. In this instance, practically, these data specify the orientation of the wall under consideration, that is its azimuthal direction. Moreover the slope of the salient also can be specified where the salients are different, for instance if the salients are natural formations. In such cases, these orientation data of the salient plane can be defined by the normal to the salient illustratively stated by the above cited 3D coordinates.

[0088] However to keep the memory module 33 compact, the orientation data of the salient can be limited to tile azimuthal data and, as regards operation, the direction of tile beam reflected at the downstream path segment is computed in light of the reflecting salients of the beam under consideration being vertical—this is the general case in an urban area.

[0089] Moreover data specifying one kind of above-ground structures, such as trees, lodgements and the like, also can be integrated into the memory module 33 of the matrix database, the computer system 30 computing the propagation conditions in the downstream path segment depending on tile specified above-ground structures.

[0090] If, in the in initial phase, the salients' edge data are integrated into the memory module 33 of the matrix database, it is feasible, in operation, to compute a direction of a refracted downstream beam on the basis of the edge data. As already indicated, refraction generally results from vertical building corners or edges. However roof edges and ridges can be specified as well in the memory module 33 in order to determine similarly a refraction direction of the incident beam showing illustratively as being deflected downward. The above deflections increase the size of the coverage zone of the cell 1 because they point the refracted beam toward a region which, in straight propagation, would be a shadow zone.

[0091] In tile initial stage, the computer system integrates attenuation data of the salients of the matrix database 33 into the memory, and, in operation, the computer system ascertains beam attenuation on the basis of the above attenuation data.

[0092] The attenuation data might relate to reflection from the salients 60 in which case the data are used to compute the attenuation of the reflected beams, which for instance might be about 7 dB, this value depending on the morphology of the building face such as glass, brick or other.

[0093] Moreover, or instead, the attenuation data might relate to the propagation through the salients 41, 42, 60, in which case the data are used to compute the attenuation of the beams which propagate in the salients as shown in FIG. 2.

[0094] In particular, the propagation-attenuation data relating to crossing the salients 41, 42 might also include data concerning a transition between propagation media and specifying attenuations of penetration into the salients or a change in the propagation medium, the data being used to ascertain a local attenuation of penetration, for instance of air/building.

[0095]FIG. 4 is similar to FIG. 3, however the building 70 shown in the former assumes, in top view, and in this instance, a simple, triangular shape for tie sake of simplification. A beam 81 is incident on a site of tie face 71 of the building 70, the site being situated in a mesh 34 completely crossed by the face 71.

[0096] According to the corresponding data of the mesh of the memory module 33 denoting the azimuthal direction of the building face 71 and furthermore indicating that the mesh 34 under consideration is completely crossed by face 71, the computer 30 determines that the incident beam 81 is reflected at the face into a beam 82 which said computer computes tile exit direction of beam 82, together with the local normal 72 to the face 71 subtends an angle which is equal and opposite to that subtended by the incident beam 81.

[0097] As indicated in this Figure, the reflected beam 82 indeed defines the main direction of a lobe 84 including accessary beams 83 included in the lobe and thereby constituting a solid angle in space: The beam 81 excites the zone on which it is incident to thereby generate a secondary source of electromagnetic radiation that is more diffuse than a conventional primary source, and has an isotropic pattern.

[0098] In order to better model the propagation of this instance, the computer system by means of an algorithm for calculating angular dispersion determines the plurality of directions of the beams 82 and 83 having specific attenuations that constitute the solid diffusion angle of the incident beam 81 beyond the site of incidence. This feature can apply to the beams reflected from inhomogeneous surfaces, such as building faces with windows and balconies, and to the diffracted beams. To illustrate the latter feature, a beam 91 is shown as being incident on tile face 71 but in a zone containing a vertical edge 73. The data of tile memory module 33 of the mesh 34 under consideration specifies the presence of the edge by indicating a main diffraction direction for an exit beam 92 and a solid angle 94 for accessary diffracted beams 93. Tile above description also applies to edges slanting with respect to the vertical. Therefore FIG. 4 also can be deemed as being a vertical section of a salient through superposed, horizontal rows of meshes, as an illustration that the diffracted beam may “crash” toward the ground within a volume that, a priori, can be assumed to be a zone of radio shade.

[0099] Be it borne in mind that the memory module 33 can simultaneously contain orientation data serving to compute a partial reflection (82) and data serving to compute a diffracted beam 92 provided that the incident beam 81, 89 has a cross-section roughly the magnitude of a mesh 34 that is only partly affected by the presence of the edge 73.

[0100] As regards FIG. 5, the roof diffraction attenuation L_(r-m), where r-m means roof to mobile (terminal), can be computed from the following formula:

L _(r-m)=−16.9−10logW+10logf+20log(hb−Hm)

[0101] or is equal to 0 if L_(r-m)<0

[0102] where

[0103] W=width of a beam received by the antenna of the base station 11

[0104] f=frequency (MHZ)

[0105] hb=height of the roof diffracting toward the mobile terminal

[0106] Hm=antenna height of the mobile terminal 21.

[0107] In order to further improve the accuracy of the estimated attenuations, the successive incidences of the beam path (51 through 55; 81 82; 91, 92) are counted in this embodiment mode and, in the event of a second incidence, in order to determine the propagation conditions beyond a second incidence point, it is assumed that the beam is polarized during the first incidence: A beam reflected or diffracted at a building face undergoes at least a partial and substantially vertical polarization. Consequently, for lack of horizontal polarization components—which were eliminated at the first incidence and which amounted to a large proportion of the total attenuation—the attenuation now has a lower magnitude at the ensuing incidences. In this way, the beam has been made to match, so to speak, the obstacles it meets. On the average, reflective attenuation changes from 7 dB for the first reflection to 3 dB at the ensuing ones.

[0108] In order to make available radio coverage data regarding the various building floors and to process the above mentioned overhang structures, the map or base terrain of attenuations is set up in three dimensions in this embodiment.

[0109] To prevent needlessly increasing computing time, tile computing system counts the consecutive beam incidences and compares their total to a high threshold value, which when reached causes operations to cease.

[0110] For the same purpose, or complementarily or instead, the computing system at each incidence determines tile cumulative attenuation and compares it to a maximum threshold value in order to discontinue operation in response to the threshold being reached.

[0111] Be it borne in mind that, in a variant of this particular embodiment tile principle of beam path reversal can be used, namely to transmit such a beam from any position of the mobile terminal 21 toward the base station 11 by shifting in the above manner the position of tile mobile terminal 21 through the entire cell 1.

[0112] However, in such a case, some uncertainty arises about the success of each test because a beam transmitted toward the base station 11 might be deflected, and inversely a beam emitted in another direction can be shifted in direction by a salient on which it is incident. Therefore a larger number of beam emitters must be used, for instance within a large solid angle containing the base station 11, and taking into account foreseeable deflections, for instance by a refracting building roof near the mobile terminal 21.

[0113] Therefore the initial position can be that intended for the base station 11 or an arbitrary position within the cells, the transmission direction of the beam is selected according to the positions of near salients and their types to cause the beam to pass near the base station 11. Accordingly a sufficient proportion of computations is available.

[0114] In particular the method of the invention can be applied to microcells in contact with or included in (macro) cells. Taking into account the substantial height of the antennas of the macrocells, propagation is hampered less by them and the computations can be carried out using a conventional propagation method. Furthermore and by means of the method of the present invention, having also computed the conditions of propagation in a microcell in contact with a macrocell, the propagation conditions in the latter are computed, followed by computational smoothing of the results of the two computations relating to a boundary zone between microcell and macrocell.

[0115] In this instance the mesh pitch of the terrain of the map 31 is about 5 m along geographic parallels and meridians. In particular it may be implemented by linear interpolation from a smaller-scale altimetry map from IGN having, meshes that are substantially large squares 50 m on each side and comprising a vector database which defines, besides the altimetry of tile terrain, the positions of the salients and what they are. The computing system meshes out said IGN map into the 5 m pitch by cutting each large square into one hundred little squares 5 m on the sides. The surface bounded by each little square in this manner determines a corresponding, sub-set of the data defining the salients positions and what they are.

[0116] Next the computing system performs smoothing or lowpass spatial filtering by using an interpolating computation that takes into account the above altimetry data from tile large squares adjacent to the one under consideration, by line, column and diagonal. Illustratively the computer 30 modulates the altitude data of the mean ground of the large square under consideration in order to ascertain therein a local, most-probable value for each little square, so as to attain in this manner a matrix sub-set of altimetric data constituting one of the plurality of the zones of tile memory module 32.

[0117] Other data of global order, for instance relating to specifying the above-ground structures such as woods, lodgements, vacant lots or other, can be computed in this manner. Data that are more specific, for instance specifying the orientation of a building face for the memory module 33, on the other hand, is preferably ascertained from surveys carried out on the terrain, for instance aerial photos. The matrix sub-set of above-ground structure data preferably, and as in this case is complemented by a matrix sub-set of the height data of above-ground structures determined from the differential of the above-ground altitude of the structure and the altitude of the ground relative to sea level. In this manner tile data of tile memory modules 32, 33 are more accurate and up to date. 

1. A method for estimating in anticipatory manner the radio coverage of a cell (1) of a cellular, wireless telephone network, by means of a radio station (11) managing tile traffic of the cell (11), using a database (32, 33) of a relief map (31) specifying said cell's positions of salients (41, 42) and what they are, said method involving a computation system dulling an operational stage, where said computation system fictitiously transmits into tile cell (1) a sampling beam that shall render the conditions of radioelectric propagation along an initial segment of a path (51) from an initial position and along one direction and in specific conditions of propagations by reading the database (32, 33) compares the path segment (51) with the data of the relief map (31) to identify tile position and the nature of a salient (41, 42) at any instantaneous incidence site of the path segment (51), determines, according the nature specified in the database of the salient (41, 42) under consideration, new propagation conditions on a downstream path segment (52) beyond said incidence site, iterates as called for the two above stages a given number of times for other downstream path segments (52-55), determines, according to tile propagation conditions along tile full path, a cumulative attenuation for any selected site of said path, and repeats all the above stages a plurality of times for a plurality of initial directions in order to sample the entire cell (1) and in this manner to ascertain a map of attenuations of the selected sites, said method being characterized in that in an initial phase, the computing system: devises tile database (32, 33) in the form of a direct-access matrix of local and independent specification of tile positions and kinds of resp. salients (41, 42) of a plurality of predetermined meshes of a corresponding mesh topology of the map (31), and stores tile geographic-orientation data of tile salients in the matrix database (32, 33), in the operational phase, the computing system: compares the position of the instantaneous site with the meshing topology to identify an incidence mesh, following incidence, computes the post-incidence propagation conditions on the basis of the local and independent salient specification data of the incidence mesh, and computes one direction of the reflected beam of a downstream segment (52) based on the orientation data.
 2. Method as claimed in claim 1, wherein storage is limited to azimuth data and, in operation, the direction of the beam reflected at tile downstream path signal (82) is calculated by assuming that tile reflection salients (71) are vertical.
 3. Method as claimed in either of claims 1 and 2, wherein data specifying the nature of the above-ground structure are integrated into the matrix database (32, 33) and wherein, during operation, the computing system computes the conditions of propagation on the downstream path segment (52) according to said nature of the above-ground structure.
 4. Method as claimed in one of claims 1 through 3, wherein, in the initial stage, salients edge data (73) are integrated into the matrix database (32, 33) and wherein, in operation, the computing system computes a direction of downstream path segment of the refracted beam (92) on the basis of said edge data.
 5. Method as claimed in one of claims 1 through 4, wherein, in the initial stage, attenuation data of the salients (41, 42) are stored in the matrix database (32, 33) and, in operation, said system computes a propagation attenuation on tile downstream path segment according to said attenuation data.
 6. Method as claimed in claim 5, wherein the attenuation data relate to reflection at the salients (41, 42) and are used to compute the attenuation on the downstream path segment (52).
 7. Method as claimed in one of claims 5 and 6, wherein the attenuation data relate to propagation through the salients (41, 42) and are used to calculate the attenuation of the beams propagating therein.
 8. Method as claimed in claim 7, wherein the propagation attenuation data through the salients (41, 42) moreover comprise data of transition between propagation media specifying penetration attenuations in the salients (41, 42) used to determine local penetration attenuation.
 9. Method as claimed in one of claims 1 through 8, wherein, with resort to an algorithm of angular dispersion computation of the beam following its incidence, said computing system computes a plurality of directions of downstream path segments having specific attenuation and constituting a solid angle (84) of beam diffusion beyond said incidence.
 10. Method as claimed in one of claims 1 through 9, wherein the computing system counts the consecutive incidences on the beam path and, in the case of a second incidence, tile said system assumes that the beam was polarized at the first incidence in order to compute said conditions of trans-propagation.
 11. Method as claimed in one of claims 1 through 10, wherein the computing system counts the consecutive incidences at the beam path and compares the total number of incidences to a high threshold value in order to discontinue the procedural stages when said threshold has been reached.
 12. Method as claimed in one fo claims 1 through 11, wherein, at each incidence, the computing system determines the cumulative propagation attenuation and compares it to a maximum attenuation threshold in order to discontinue the procedural stages when said threshold has been reached.
 13. Method as claimed in one of claims 1 through 12, wherein tile attenuation map is stored in 3D.
 14. Method as claimed in one of claims 1 through 13, wherein, to devise tile database (32, 33), tile computing system represents the map (31) by a bundle of vertical-extension pixel links and selects pixel links to constitute stacked volume elements of volume meshes each comprising particular data.
 15. Method as claimed in one fo claims 1 through 14, wherein the computing system assumes a position assigned to the base station (11) being the initial position.
 16. Method as claimed in one of claims 1 through 14, wherein the computing system assumes that the initial position is an arbitrary position within the cell and wherein the beam transmission direction is selected according to the positions and the kinds of structures in the vicinity in order that said beam shall travel near the base station (11).
 17. Method as claimed in one of claims 1 through 16, wherein, after having computed in this manner the conditions of propagation within a microcell in contact with the cell, the computing system calculates propagation conditions within the cell and next computes smoothing the results of two computations relating to a boundary zone between the said cell and microcell. 